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ABSTRACT 

In recent years, near-IR and X-ray flares have been detected from the Galaxy's central radio 
point source, Sagittarius A* (Sgr A*), believed to be a ~ 3 x 10 6 M Q supermassive black hole. 
In some cases, the transient emission appears to be modulated with a (quasi-)periodic oscillation 
(QPO) of <~ 17 — 20 minutes. The implied ~ 3r$ size of the emitter (where rg = 2GM/c 2 
is the Schwarzschild radius) points to an instability — possibly induced by accretion — near the 
Marginally Stable Orbit (MSO) of a slowly spinning object. But Sgr A* is not accreting via a 
large, 'standard' disk; instead, the low density environment surrounding it apparently feeds the 
black hole with low angular momentum clumps of plasma that circularize within ~ 10 — 300 r$ 
and merge onto a compact, hot disk. In this Letter, we follow the evolution of the disk following 
such an event, and show that a Rossby wave instability, particularly in its magnetohydrodynamic 
(MHD) form, grows rapidly and produces a period of enhanced accretion lasting several hours. 
Both the amplitude of this response, and its duration, match the observed flare characteristics 
rather well. 
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1. Introduction 

The prospect of identifying phenomena within 
several Schwarzschild radii (rs = 2GM/c 2 ) of Sgr 
A* has improved considerably in recent years, with 
the detection of near-IR (Genzel et al. 2003) and 
X-ray (Belanger et al. 2005) flares modulated with 
an average period of s=s 17 — 20 minutes (Belanger 
et al. 2006). This (quasi-)period corresponds to a 
modulation near the MSO of the black hole (Melia 
et al. 2001; Liu & Melia 2002), possibly suggest- 
ing a transient event associated with a wave pat- 
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tern co-rotating with the gas, or perhaps a "hot" 
spot where matter has fallen in from larger radii, 
impacting the small disk (Falcke & Melia 1997). 
In this Letter, we examine numerically the Rossby 
wave instability (RWI) induced by a clump of mag- 
netized plasma merging onto the disk and focus on 
its implications for the subsequent flow of matter 
through the MSO toward the event horizon. 

For a black hole mass M ~ 3 x 10 6 M , the 
radius r w 3r$ inferred from the average period 
is interesting for several reasons. Aside from con- 
straining the flaring event to a region very near the 
MSO, it supports the view that there is indeed a 



hot, magnetized disk confined to the inner ~ lOrg 
that can simultaneously produce a ~ 10% polar- 
ized component and a position angle that rotates 
by 90° across the mm/sub-mm bump (Bromley, 
Mclia, & Liu 2001). In addition, strong variability 
is seen near the middle of these events, in which 
the X-ray flux drops by a factor of 40-50 % in 10- 
15 minutes. Simple light travel time arguments 
constrain the emitting region to be no bigger than 
- 17 - 34 r s (Melia 2001). This tight orbit sug- 
gests that the X-ray and mm/sub-mm photons 
are produced by the same medium, probably via 
synchrotron-sclf-Compton processes. Finally, this 
is in line with the measured intrinsic size of Sgr 
A* at A7 mm (Bower ct al. 2004). 

Extensive hydrodynamic (Melia & Coker 1999; 
Cuadra et al. 2005) and MHD (Igumenshchev 
& Narayan 2002) simulations show that for the 
stellar-wind fed conditions at the Galactic center, 
the average specific angular momentum of gas cap- 
tured gravitationally by Sgr A* is too small to sus- 
tain a 'conventional' (i.e., typically large ~ 10 5 r$) 
disk. Instead, only clumps of plasma with rela- 
tively small angular momentum venture inwards 
and merge with — essentially, 'rain' onto — the com- 
pact disk at the circularization radius, which for 
Sgr A* is < 10 - 10 3 r s (see also Melia & Falcke 
2001, Melia 2006). In this paper, we model the 
transient disruption to the disk induced by blobs of 
plasma circularizing within the inner ~ 10 — 20 rs, 
and show that the ensuing dynamical evolution 
has characteristics in common with those inferred 
from the flares. 

2. Magnetic Rossby Wave Instability 

The instability we are here concerned with has 
a long history, dating back to Lovelace & Hohlfcld 
(1978), who showed that a disk presenting an ex- 
trcmum of a quantity £ (later dubbed vortensity) 
was subject to a local instability of Rossby vor- 
tices. The requirement of an extremum is similar 
to that giving rise to the Kclvin-Hclmholtz insta- 
bility of sheared flows. More recently, Lovelace 
ct al. (1999) (see also Li ct al. 2000, Li et al. 2001) 
renamed it the Rossby Wave Instability (RWI) and 
developed the theory and numerical simulation. 
Varniere & Tagger (2005) used it to address the 
presence of a 'Dead Zone' in protostellar disks. 

In isothermal, unmagnetizcd disks, £ is the spe- 



cific vorticity averaged across the disk thickness, 

S 2fi£ ' [ > 

where E is the disk's surface density, ft its rota- 
tion frequency, and k 2 = 4ft 2 + 2QQ'r is the 
epicyclic frequency squared. The extremum of £ 
appears to be due to an extremum in the radial 
density profile. However, starting from a very dif- 
ferent point of view — the seismology of relativis- 
tic accretion disks — Nowak & Wagoner (1991) de- 
scribed g-modes that intoduce an extremum in k, 
and thus also in £ (though this time due to the 
relativistic rotation curve), near the MSO (where 
k vanishes). The g-mode of Nowak & Wagoner is 
thus similar to the RWI, but they seem to have 
missed the mechanism that makes it unstable. Its 
tentative association with a QPO in X-ray binaries 
does not seem to be supported by observations. 
We will see, however, that it plays an important 
role in what follows. 

Many of the details regarding how the instabil- 
ity is driven have already appeared in the work 
cited above. The main idea is that Rossby waves 
in disks form normal modes trapped near the ex- 
tremum of £. The waves propagate either outside 
or inside their corotation radius, i.e., the radius 
where the wave rotates at the same frequency as 
the gas, depending on the sign of the gradient in 
£ (Tagger 2001). They have a positive energy in 
the first case, and a negative energy in the sec- 
ond. In a forthcoming paper, we shall discuss 
the MHD form of the RWI, for a disk threaded 
by a vertical (poloidal) magnetic field B (r). Its 
properties are essentially the same as those dis- 
cussed above, except that here the critical quan- 
tity is £b = K 2 E/(2ri_Bg), and the growth rate 
can be higher because of the long-range action 
of the Lorentz force coupling the Rossby vortices 
(Tagger & Pellat 1999). Another difference is that 
the energy and angular momentum of the Rossby 
waves can escape the disk vertically as Alfven 
waves propagating along the magnetic field lines 
(Varniere & Tagger 2002). 

Let us finally mention that even without an 
extremum in £ or £g, normal modes may still 
occur if the disk has an inner edge that reflects 
spiral density waves. In magnetized disks this 
gives the Accretion-Ejection Instability (Tagger 
& Pellat 1999). Recently, Li et al. (2003) dis- 
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Fig. 1. — Surface density (top) and (bottom) 
profiles, at startup (full) and after 20 (dashed) and 
250 (dot-dashed) MSO orbital times in our first 
simulation. Radii are scaled to tmso- 

cussed the non-axisymmetric g-modes and con- 
cluded that those without nodes in their vertical 
structure (that we analyze here) are not unstable. 
But this conclusion, based on the absence of a re- 
flecting inner edge in a relativistic disk, neglected 
the possibility of trapped modes based on Rossby, 
rather than density, waves. 

3. Numerical simulations 

The simulations presented here are relevant to 
the conditions in Sgr A*. We use the 2-D MHD 
code introduced in Caunt & Tagger (2001), that 
considers an infinitely thin disk in vacuum. This 
is well adapted to the physics we consider, since 
the RWI and spiral modes have a structure that 
is essentially constant across the disk thickness, 



allowing us to integrate the dynamical equations 
over z. (On the other hand, this simplification 
forbids us to study the Magneto- Rotational Insta- 
bility (MR1), whose structure does vary with z.) 

The radial grid points arc distributed loga- 
rithmically, allowing simulations that extend far 
in radius (typically with r max /r min « 50), so 
as to minimize unwanted boundary effects. We 
have modified the code by including the pseudo- 
Newtonian potential of Paczinsky-Witta, $(r) = 
GM/(r — rs)- This allows for the existence of an 
MSO, tmso — 3rs, where k vanishes. Further- 
more, at the inner grid point the gas is allowed 
to cross freely during its final collapse toward the 
central object. The simulations are carried out on 
a polar grid with n r — 256 and ng = 128 points, 
and extend from ~ 0.6 to 30 tmso- The MSO is 
at the 30th radial grid point, allowing us to resolve 
properly the gas flow (including the sonic point) 
in the plunging region; the infall velocity at the 
inner boundary is typically 3 to 5 times the sound 
velocity. Since pressure does not play a major role 
in the physics of Rossby waves, we fix the sound 
speed as a function of radius, which for simplicity 
is kept constant throughout the simulations. We 
first present the simulation of a disk with initially 
a power-law surface density profile, £ ~ r -1 / 2 , 
and a sound speed eg = .05 rQ(r), (giving a con- 
stant disk aspect ratio h/r). The initial magnetic 
field is such that (3 = 8irp/B 2 = 2 throughout the 
disk. After a few orbits, the density profile near 
the MSO reaches a quasi-steady state, as shown 
in Figure 1 (here and below we scale all radii to 
I'm so and times to the orbital period there). This 
is very similar to what is obtained with 3-D codes 
(Hawley & Krolik 2001; Machida & Matsumoto 
2003). An m=l perturbation, identified as a g- 
modc instability by the RWI mechanism, sets in, 
gradually eroding the inner edge of the disk (since 
there is no turbulent accretion from farther out). 
This builds up the extremum of Cb just beyond 
the MSO. The instability has its corotation radius 
near the extremum, as expected. The accretion 
rate M through the disk peaks when the instabil- 
ity gets established, and then decreases gradually 
for the rest of the simulation, as the density max- 
imum slowly moves outward. 

In our second calculation, we simulate a flare 
event by adding to the 'standard' disk a strong 
bump in density and magnetic field centered at 
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4 tmsOi where it generates a minimum in 
This bump represents 56% of the total mass in 
the central 6 tmso ■ Pursuant to our discussion in 
the introduction, we assume that a blob of magne- 
tized plasma with very little angular momentum 
has merged onto the disk, rapidly circularizing to 
form a ring at this radius. As shown in Figure 2, 
the evolution is much more violent. In a nearly 
sudden episode between 40 and 60 MSO orbital 
times (hardly 2 periods at ^mso, or in real time 
about 5 hours), most of the added gas has been 
transported to the inner region of the disk by vi- 
olent instabilities developping over the whole re- 
gion. 

We see m = 1 modes developing (as previously) 
near the inner edge, and m = 2 in the bump re- 
gion, with a corotation near the extremum of £g. 
These modes appear to be non-linearly coupled; 
similar coupling has been observed in simulations 
and observations of spiral galaxies (Tagger ct al. 
1987; Masset & Tagger 1997). This coupling can 
be explosive (i.e., result in faster than exponential 
growth) when it involves modes of positive and 
negative energy (Dikasov et al. 1969). Figure 3 
shows the surface density profile at t = 40£mso> 
when the density bump starts being disrupted by 
the instability. Figure 2 shows that after this rapid 
phase, the radial profile of becomes perfectly 
flat beyond the MSO region, confirming qualita- 
tively that it is that profile to which the instability 
reacts. 

Producing realistic light curves involves much 
more physics than we include here, and is be- 
yond the scope of the present simulations. How- 
ever, figure 4 shows M through the inner edge of 
the simulation box. Notice the brief episode of 
strong, enhanced accretion between t = 40^mso 
and t = GOt mso, with the largest amplitude oc- 
curring over a time interval of roughly Hmso (i-e., 
~ 140 minutes). The similarity between this time 
and the duration of a typical flare (~ 2—3 hours) is 
very encouraging. Also, there is clear evidence of 4 
or 5 quasi-periodic pulsations during this phase of 
enhanced accretion, mimicking the behavior wit- 
nessed during several IR flares and the XMM event 
of 2004. This appears to be due to an axisym- 
metric g-mode, trapped near the extremum of £g 
and non-linearly excited by the instabilities. The 
details may vary with the numerical setup, e.g., 
with the amount of artificial viscosity; this is rea- 




Fig. 2. — Surface density (top) and £j$ (bottom) 
profiles, at startup (full) and after 20 , 40, 47, 60, 
75 MSO orbital times during our second simula- 
tion. 



sonable since this is the only dissipation limiting 
the growth of the instability in the simulations 
(besides the destruction of the bump itself), but 
the main features (duration of the flare and pul- 
sations) remain the same. 

Simulations of less magnetized or unmagnetized 
disks, or of a weaker bump, give a very similar, 
though somewhat slower, evolution. Placing the 
bump at 6 or 8 tmso delays the development of 
the instability, but does not change the duration 
or the amplitude of the 'flare', which occurs only 
when the mass of the accreted gas approaches the 
MSO. However this might result in changes of the 
mm and cm cmissivity relative to the IR/X-ray 
emission, so a comparative study of the broadband 
spectrum of many flares may provide indirect ev- 
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Fig. 4. — Accretion rate through the inner edge of 
the simulation box in the second simulation. 



Fig. 3. — Surface density in our second simulation 
at t = 40 tMSO, when the bump at 4 tmso (red, 
dark) starts being disrupted by the instability. 

idence of the variation in clumping radius. 

4. Conclusion 

The RWI, particularly in its MHD form, can 
produce the behavior required to account for the 
near-IR and X-ray flares in Sgr A*. Unstable both 
near the MSO and in the annulus where a clump of 
interstellar plasma has merged onto the disk, it de- 
velops as strong Rossby vortices whose non-linear 
evolution can be remarkably fast. The instabil- 
ity grows in a few rotation times, and causes the 
clump to accrete to the inner edge of the disk. The 
resulting period of enhanced M at the inner edge 
corresponds to the final stage of this collapse, so 
that its duration (comparable with the observed 
duration of the flares) does not depend on the 
initial location of the clump. Furthermore, the 
enhanced M is accompanied by QPO very sim- 
ilar to those observed both in IR and in X-rays. 
The quasi-period reflects the underlying Keplerian 
frequency near the MSO, though the exact quasi- 
period may be dependent on our use of a pseudo- 
Newtonian potential, rather than a fully relativis- 
tic model. Natural extensions will thus be to in- 
clude the metric of a spinning Black Hole, and to 
consider the relevance of the RWI to explain the 
high-frequency QPO in microquasars. 



While current wisdom seems to be settling on 
an accretion instability for the origin of Sgr A*'s 
flares, the process by which the actual emission 
occurs is uncertain. The favored scenario right 
now connecting the mm/sub- mm to NIR and X- 
rays portions of the spectrum (Eisenhauer et al. 
2005; Liu & Melia 2001) is that the former is due to 
synchrotron, whereas the X-rays are produced by 
synchrotron-self-Compton (SSC). It turns out that 
producing the right blend of physical conditions to 
fit both the NIR and X-ray flare emission (under 
the assumption that the two occur more or less 
simultaneously) is not trivial. 

In a related paper (Liu, Melia & Petrosian 
2005), we considered a scenario in which the radia- 
tive manifestation of the flare occurs when the en- 
ergy released by the accretion instability is trans- 
ferred to a population of nonthermal particles via 
stochastic acceleration by the turbulent compo- 
nent of the magnetic field. Indeed, as mentioned 
in §2, Rossby waves in the disk extend into the 
corona as Alfven waves, carrying much of the 
accretion energy. Clearly, much work remains 
to be done in coupling these complementary ap- 
proaches. One of the main results reported in 
this paper is that the enhanced M during a flare 
event is consistent — both in amplitude and tem- 
poral evolution — with the characteristics observed 
in Sgr A*'s flares. We are in the process of self- 
consistently calculating the particle acceleration 
following an RWI-induced event, and will report 
the results of this work in a future publication. 
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